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Quantum mechanical entanglement is a resource for quantum computation, quantum teleporta- 
tion, and quantum cryptography. The ability to quantify this resource correctly has thus become 
of great interest to those working in the field of quantum information theory. In this paper, we 
show that all existing entanglement measures but one fail important tests of fitness when applied 
to n particle, m site states of indistinguishable particles, where n, m > 2. The accepted method 
of measuring the entanglement of a bipartite system of distinguishable particles is to use the von 
Neumann entropy of the reduced density matrix of one half of the system. We show that expressing 
the full density matrix using a site-spin occupation number basis, and reducing with respect to that 
basis, gives an entanglement which meets all currently known fitness criteria for systems composed 
of either distinguishable or indistinguishable particles. 

We consider an output state from a previously published thought experiment, a state which is 
entangled in both spin and spatial degrees of freedom, and show that the site entropy measure gives 
the correct total entanglement. We also show how the spin-space entanglement transfer occurring 
within the apparatus can be understood in terms of the transfer of probability from single-occupancy 
to double-occupancy sectors of the density matrix. 

PACS numbers: 03.67.-a,03.65.Ud,05.30.-d 



I. INTRODUCTION 



The peculiarly non-local correlations exhibited by the 
states of quantum systems are key to the implementation 
of quantum information processing technologies, such as 
quantum computation and quantum teleportation. How- 
ever, it is easily shown that the correlations due to the 
(anti)symmctrization of the states of indistinguishable 
bosons (fermions) are not themselves a physically use- 
ful resource for quantum information technologies: for 
example, there is no measurement we can make locally 
on a fermion in a localized state which is affected by 
the existence of identical fermions in other parts of the 
universe jlj|. However, it is possible to produce entangle- 
ment that is a resource for QIT by suitable preparation: 
for example, by producing a |^)~ Bell state of the spins 
of two fermions. Indeed, in practice, many potential im- 
plementations of QIT involve identical particles (such as 
photons, electrons, or protons) as 'carriers' of entangle- 
ment. It is therefore important to be able to quantify the 
degree of 'useful' entanglement in a system of identical 
particles. 

Discussion of the entanglement between pure states 
of indistinguishable particles has previously been dealt 
with almost as a separate topic from that of distinguish- 
able particles. It is the aim of this paper to show that 



the entanglement of pure states of either type of particle 
can be described within the same theoretical framework. 
This framework involves the von Neumann entropy of 
the reduced density matrix for the subsystem whose en- 
tanglement with the rest of the system we wish to find, 
expressed in an occupation number basis Q. This also 
allows us to understand better the division between spin 
and spatial entanglement in systems where both may ex- 
ist, and the manner in which entanglement may be trans- 
ferred between spin and space. It is important to em- 
phasize that we consider in this paper only pure states of 
the full system. It is already known that for such states 
the von Neumann entropy provides the correct measure 
of entanglement between two distinguishable subsystems 
||. We do not address the case of an overall mixed state, 
for which the definition of an entanglement measure is 
more subtle 101. 
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In section II we discuss the partitioning of Hilbert 
space that is implicit to any meaningful definition of en- 
tanglement. In section [II we review some requirements 
for a successful entanglement measure, and consider the 
extent to which three potential definitions meet these re- 
quirements. In section ^ we show that Zanardi's site 
entropy measure passes all the tests, and can be related 
to the conventional definition of entanglement in the limit 
where the exchange symmetry of the particles is irrele- 
vant. Finally, in section [y] we use Zanardi's measure to 
discuss spin-space entanglement transfer. 
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II. METHODS OF PARTITIONING HILBERT 
SPACE OF TWO ENTANGLED SPINFUL 
PARTICLES 

Implicit to any measure which attempts to describe the 
entanglement of two subsystems is an assumption about 
the correct manner in which to partition the total Hilbert 
space. In this section we consider the requirements for 
a correct partitioning, and look at how this is actually 
performed by existing entanglement measures. We will 
frequently need to talk about the states of internal de- 
grees of freedom of particles. Therefore, for brevity we 
will henceforth refer to any states of such internal degrees 
of freedom simply as 'spin states'. 



1. Requirements for partitioning 

Tensor product structure In order to express entan- 
glement between two components of an entangled system, 
some kind of partitioning of their Hilbert space is neces- 
sary in order to identify the 'components'. Our aim is to 
quantify the entanglement resource shared between parts 
A and B of a composite quantum system. These parts 
may be identified with particles (in the case of a state of 
the system where the particles are localized), with sites 
(in the case of a state of the system where the particles 
are delocalized over sites) , or with some arbitrary subdi- 
vision of an experimental apparatus (as we will examine 
in section VA). For the purposes of the greater part 
of this paper, we will consider subsystems of a system 
as being synomymous with sites. But it is important to 
emphasize that our conclusions are more general: they 
apply to any division of a system into subsystems. 

For entangled states of distinguishable particles (or 
particles which are effectively distinguishable because of 
their localization) we would normally use a tensor prod- 
uct structure H — Ha <8> Hb where Ha and Hb are 
Hilbert spaces for states of particles in parts A and B. 
It is important that we correctly partition the Hilbert 
space because this ensures that basic operations such as 
the partial trace pb = trsp are valid. The partial trace 
is the correct and only way to describe the properties of 
one part of a composite quantum system when nothing is 
known about the other parts, as it gives the correct mea- 
surement statistics for observations on that subsystem 



But if we try to use the tensor product structure parti- 
tioning for entangled states of indistinguishable particles, 
we run into two problems: 

• The Hilbert space of two indistinguishable parti- 
cles is a symmetric or antisymmetric product, not 
a direct product. 

• There is no correspondence between the particles 
and the subsystems used in the partitioning. 



Delocalization For spin-only entangled states of dis- 
tinguishable particles — i.e. states where we have unam- 
bigously given one particle to Alice, and the other to 
Bob — the phrase 'the states of Alice's spin' is completely 
equivalent to the phrase 'the states of Alice's particle'. 
There is no ambiguity about which particle Alice has in 
her possession at any time, and therefore there is no logi- 
cal difference between a one-site (local) unitary transfor- 
mation, and a one-particle (possibly non-local) unitary 
transformation. Thus when deciding on a basis in which 
to describe the spin-only entanglement of a system of 
distinguishable particles it may seem a matter of taste 
whether spin states should be assigned to particles, or to 
sites. 

However, it is perfectly possible to write down states in 
which each particle is shared between Alice and Bob. An 
example of such a 'spin-space entangled state' is obtained 
if we put particle 1 into -i=(A | +B |) and particle 2 into 

-^{A i +B I), where A, B are site labels. 

Indistinguishability When the entangled particles are 
indistinguishable, we can no longer be sure which particle 
Alice has in her possession. The distinction between one- 
particle unitary transformations, and one-site unitary 
transformations becomes relevant. Entanglement should 
be invariant under one-site unitary transformations, but 
not necessarily under one-particle unitary transforma- 
tions, which may generate entanglement if they involve 
both subsystems. An entanglement measure which works 
successfully for indistinguishable particles must respect 
this distinction. 

The natural way to achieve this distinction is to use 
a basis which assigns spin states to sites rather than 
particles. 



2. Partitioning used by existing entanglement measures 

When partitioning the total Hilbert space of two en- 
tangled quantum systems, we need to ask ourselves: 

• For indistinguishable subsystems: to what ex- 
tent can my system be regarded as a symmet- 
ric/antisymmetric product of the single-subsystem 
states? 

• For distinguishable subsystems: to what extent can 
my system be regarded as a direct product of the 
single-subsystem states? 

In most descriptions of entanglement, the tensor prod- 
uct structure is used, for example in the entanglement 
measure introduced by Wootters ||. This measure is 
suitable for spin-only entanglement of localized distin- 
guishable particles. However, it does not describe which 
site a particle occupies, so is not suited to describing 
either entangled indistinguishable particles, or entangled 
states of distinguishable particles where the 'particle' and 
'subsystem' divisions do not coincide. 
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One example where indistinguishable particles have 
been treated is by Schliemann et al. @ [§ , who explicitly 
consider the antisymmetric product space belonging to 
two fermions, each of which inhabits a four-dimensional 
one-particle space. They write a general state in the six- 
dimensional two-particle Hilbert space as 



\w) 



a,be{l,2,3,4} 



w abC' a C b 



Ulo> 



(1) 



where a,b run over the orthonormahzed single particle 
states, and Pauli exclusion requires that the 4x4 coeffi- 
cient matrix w is antisymmetric: w a b = —Wb a - 

It may seem that Schliemann's partitioning is indeed 
in terms of sites rather than particles, since the single 
particle states are labelled by sites. But, as we shall 
see later, Schliemann's measure is derived by considering 
the number of elemental Slater determinants needed to 
expand the entangled state. It is therefore actually a 
particle-based, rather than a site-based, description of 
entanglement. As a consequence, as will be shown later 
in this paper, it suffers from a number of serious flaws; 
in particular, it is possible to devise one-site (i.e. local) 
transformations which generate 'entanglement' according 
to the Schliemann measure. 



III. REVIEW OF EXISTING ENTANGLEMENT 
MEASURES 

A. Desirable properties of any entanglement 
measure 

What are the desirable properties of an entanglement 
measure? 

Invariance under local unitary transformations. If a 
measure is correct, it should not be possible to generate 
'entanglement' using only unitary transformations local 
to a particular site. 

Non-invariance under non-local unitary transforma- 
tions. Conversely, it should be possible to find non-local 
(i.e. multisite) unitary transformations which change the 
entanglement. 

Correct behaviour as distinguishability of subsystems A 
and B is lost. A correct measure should reflect the fact 
that entanglement is affected when the distinguishability 
of the subsystems involved is lost. A simple example of 
this is as follows. For two fermions whose spin degrees of 
freedom are maximally entangled, we require that as the 
overlap of the single-particle spatial wavefunctions ap- 
proaches unity the entanglement should asymptotically 
approach zero. This is easily seen by considering the full 
expression for the Bell basis states in terms of Slater de- 
terminants. 

If the two fermions are localized, one in site A and one 
in site B, then the l^)* Bell state can be written 



|*) = 



1 



V2 



(IU)±UT» 



(2) 



where the full expression for ||J.) is 



IU) 



1 



<m(i)T(i) Mi)i(i) 

4a (2) T (2) <M2) I (2) 



(3) 



When the two fermions are brought together to occupy 
the same site, the spatial parts of the two single-particle 
states coincide, i.e. <f>A — > 4>b, and we have 



IU) - -m 



l 

V2 



0(1) t (1) 0(1) 4 (1) 
0(2) t (2) 0(2)4(2) 

0(1)0(2) 



T(l) 
T(2) 



(1) 
(2) 



(4) 



where is the same spatial state for sites A and B. 

A similar result is obtained for |J,|), but with an ex- 
change of columns and therefore the same result applies 
for it as for \\\) but with an overall minus sign. Hence 
|44) and | J, 4) are now linearly dependent, and the be- 
haviour of the entangled Bell state is: 



I* 



± <t>A— ><t>E 



1 0(1)0(2) 



V2 

( 



^/2 

r(i) 

T(2) 4(2) 



± 



1(1) 

1(2) 



T(i) 

4(2) 



i 0(1)0(2) ( 4(1) 4(2)- 4(1)4(2) 



±(4(1)4(2)- 4(1)4(2)) 



and hence 



I*) 



0(1)0(2) 



T(i) 1(1) 

T(2) 4(2) 



(5) 



(6) 



up to a normalization factor, whereas |\&) + — * because 
of Pauli exclusion. 

Thus the one ebit of entanglement present in a \&) 
state should be destroyed as the spatial overlap of the 
two fermions' wavefunctions asymptotically approaches 
unity — in the case of | V E') + because the state itself is de- 
stroyed, and in the case of \^)~ because the entangled 
Bell state becomes a non-entangled product state. (At 
least, this is the case if neither Alice nor Bob can mea- 
sure with spatial resolution sufficient to determine the 
substructure of the spatial state 0.) A correct entangle- 
ment measure should reflect this fact. 

For a pair of bosons in the |vE') ± state, exactly the 
same loss of entanglement would occur, although the be- 
haviours of and l^) - are exchanged, due to the 
change of sign introduced by the use of permanents rather 
than determinants. 
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B. Wootters measure for distinguishable particles 
(tangle) 

Wootters ||, ||, considers a particular state of two 
distinguishable particles: 



= o|TT> + 6|U) 

|a| 2 + |6| 2 



c|4T) 

l„|2 , 



dill), 



Ml 2 



1 



(7) 



where it is implicit that each particle occupies a definite 
spatial state. Since the tensor product decomposition 
of H allows us to define a reduced density matrix ps 
describing the mixed state of system B, the von Neumann 
entropy of ps is a natural measure of entanglement. The 
Wootters entanglement is simply a reexpression of the 
von Neumann entropy of pb, and is defined as 



.1 



where 



E = h[-(i + VT^P)}, 



h(x) = -(xlog 2 x + (1 - aj)log 2 (l - x)) 



(8) 



(9) 



and the quantity r is known as the 'tangle' and is defined 
by 



= 4\ad-bc\ 2 



(10) 



Since -E-woottcrs expresses the entropy of a single site, 
there is no single-site operation which can affect it. 

The Wootters measure applies only when the particles 
are totally distinguishable by virtue of occupying distinct 
sites. But our aim is to describe more general states in 
which each particle occupies a superposition of sites — 
what happens if we simply go ahead and use the Wootters 
measure regardless? Since the Wootters measure does 
not depend on the nature of the spatial states, there is 
no way its value can change. So for example, there is 
no way that the Wootters entanglement of a | V I') ± Bell 
state will ever be affected by the spatial overlap of the 
single-particle wavefunctions of the constituent particles. 



C. Schliemann measure for fermions 

Schlicmann et al. Q define the entanglement of spin 
states of a pair of fermions by 



T)(w) := \(w\w)\, 
where the dual w of w is defined by 

w ab = \e abcd w cd 
and the inner product is expressed as 

^ w* 6 Wcd(0|cbC o 44|0) 

abed 
abed 



(ii) 



(12) 



(w\w) 



= 8(wi 2 U>34 + W13W42 + W14W23)- (13) 



A similar definition was introduced for a pair of bosons 
by Paskauskas and You |0. 

Slater decomposition form. It is possible to relate the 
Schliemann measure n to the number of elementary Slater 
determinants that are required to construct the entangled 
state. The Hilbert space for a two-fermion, K-site system 
is the antisymmetric space A(C 2K ®C 2K ). Any vector in 
this space can be represented in terms of single particle 
functions f a u\ > which are members of the single-particle 
space C 2K , by the Slater decomposition 



I*) 



1 



K 



K 



Zi/Iim/LjO) (14) 



The number of non-zero coefficients required to con- 
struct \^f), i.e. the number of elementary Slater deter- 
minants, is known as the Slater rank of the entangled 
state. Then for a two-fermion two-site system, |\&) has 
Slater rank 1 (consists of a single Slater determinant) iff 
»7(l*»=0. 

Behaviour as overlap of single particle wavefunctions 
is increased. This entanglement measure behaves cor- 
rectly as the overlap is increased between the single- 
particle wavefunctions of the particles, as is shown in 
Figure [l} 

Relation to Wootters measure. Let us consider how 
the Schliemann measure works for the class of states con- 
sidered by Wootters: 



\<f>) = a\n) 



b\n) 
+ \b\ 2 



c|lT>+d|U>, 

Ic| 2 + |rf| 2 = l 



(15) 



In the representation used by Schliemann, we can write 
the w-matrices for the two-particle basis states in the 
A\,Al,B1,Bl basis as 
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(16) 



Therefore the state considered by Wootters, 

|0) =a|TT)+6|U> + c|iT> + d|U>, 

|a| 2 + |6| 2 + |c| 2 + \d\ 2 = 1, has the coefficient matrix 




(17) 



and thus we obtain the relation 
\(w\w)\ 



\e abcd w ab w cd \ 



|8(u;i2l«34 + + Wl4W>23)| 

2\ad- bc\ = sfr. 



(18) 
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Hence for the state of two distinguishable particles con- 
sidered by Wootters, the Schliemann measure rj is related 
to the Wootters 'tangle' r by 



V 



(19) 



Non-invariance under local unitary transformations. 
We can however easily show that there are local (one- 
site) unitary transformations that generate 'entangle- 
ment' by the Schliemann measure. Consider this two 
particle state: 

m = ^(4 T + 4 T )^(4 + c k)io>. (20) 

The physical interpretation of this state is that it de- 
scribes a doubly filled 'molecular orbital' 



\A) + \B) 
V2 



(21) 



where \A), \B) are the spatial states for sites A, B re- 
spectively. 

Its antisymmetric coefficient matrix is 



(22) 



giving a Schliemann entanglement of r\ — (no entangle- 
ment, since it is a single Slater determinant). 

Now consider applying the infinitesimal one-site two- 
particle unitary transformation (1 — ieH) with H = 
UnA^riAi- This purely local operation transforms the 
ui-matrix to 





(23) 



which gives a Schliemann entanglement of r) = 8| ~" | = 
2eU which is non-zero to first order in e. We have suc- 
ceeded in generating Schliemann 'entanglement' via a 
purely local unitary operation, something that it should 
not be possible to achieve. 

Invariance under non-local unitary transformations. 
Now consider an infinitesimal two-site one-particle uni- 
tary transformation. We would expect such a transfor- 
mation to lead to a change of entanglement, yet we can 
construct an example under which the Schliemann mea- 
sure is invariant. Our example is generated by a Hamilto- 
nian describing intersite hopping accompanied by a spin- 
flip: 



H = t(c At c Bi + cLc At ). 



(24) 



(The spin-flip is introduced so that our state is not an 
eigenvector of H). The Hamiltonian's action on our ex- 
ample state is 



H\$) 



(25) 



Hence applying the infinitesimal unitary transformation 
(1 — ieH) with this operator to our example state \ip) = 
'~S( c Ill + c bi)|0)j we obtain a w matrix 



with extra terms ±iet in the A\ B] and B[ A{ locations: 



/ 

-1 

-iet 

V-i 



1 


—iet 


1 





-1 


—iet 


1 







iet 


-1 






(26) 



which has a Schliemann entanglement 



V = = 0(e 2 )- 



(27) 



Thus, to first order in e, the Schliemann entanglement 
of our example state is unaffected: even though the trans- 
formation introduces new correlations between the spin 
states of the particles on sites A and B. 

Understanding the anomalous behaviour of the Schlie- 
mann measure in terms of the Slater decomposition. 
The Slater decomposition representation of an entan- 
gled two-fermion, two-site state described earlier pro- 
vides a particularly simple way of understanding why the 
Schliemann measure does not behave correctly under ei- 
ther two-site one-particle or one-site two-particle unitary 
transformations. 

According to Schliemann et al, a two-fermion two-site 
state is entangled iff it has a Slater rank greater than one. 
It is well-known that a one-particle unitary transforma- 
tion applied to a Slater determinant will produce another 
Slater determinant, whereas a two-particle transforma- 
tion will produce a superposition of Slater determinants. 
Therefore any one-particle two-site unitary transforma- 
tion will not affect the Slater rank of a state and so will 
not change the Schliemann entanglement, despite being a 
non-local transformation. Similarly, all two-particle one- 
site unitary transformations will modify the Slater rank 
of a two-fermion two-site state, and therefore will change 
the Schliemann entanglement, even though they arc lo- 
cal. Schliemann's measure therefore fails to behave as 
we expect. The entanglement measures introduced in 
ifTTf and |Q suffer from analogous problems, since both 
are based on the rank of the state. 



D. Zanardi measure 

Zanardi [g] considers the Fock space of N spinless 
fermions in a lattice with L sites. The state space Hl(N) 
for this system is given by 



H L {N) :=span{|A)/AGPi v } 



(28) 



where the antisymmctrized state vector \A) is given by 
the Slater determinant 



1^4) 



1 



(-1) 



I -PI 



(29) 
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and where Pr denotes the family of N-site subsets of the 



site labels 
„-th 



L 

and \ipj Pm ) is the single particle state for the 
j site where j is a member of the subset . 

For some G Hl{N), the local density matrix for 
the j th site is given by 



giving 



Pi :=tril*)<*| 



(30) 



where try denotes the trace over all but the j site, and 
therefore the von Neumann entropy of pj is a measure of 
the entanglement of the j th site with the remaining N-l 
sites. We will now show that, unlike the other candidates, 
Zanardi's measure possesses all the desirable features of 
an entanglement measure that we have listed above. 



IV. 'SITE ENTROPY' ENTANGLEMENT 
MEASURE 

A. Application to an example state 

Let us now investigate further the properties of Za- 
nardi's 'site entropy' entanglement measure. Our conclu- 
sions are that using a binary site-spin occupation num- 
ber basis for the full density operator for an entangled 
system, and then reducing the density operator with re- 
spect to this basis, gives a reduced density matrix whose 
von Neumann entropy appears to be a correct measure 
of entanglement under all circumstances, and for all spin 
statistics. This is due to the fact that Fock space (to 
which this representation maps the Hilbert space of a 
set of indistinguishable particles) has a natural product 
structure. 

For example, for the (fermionic or bosonic) state con- 
sidered in a previous section 

iv;} = ^(c% + 4 T )^(4 + 4 x )io), (3i) 

the density operator for the full system is 

P= T( C AT + C B])( c M + c Bl)\°)(°\( c Al +CBl)(c A ] +C BT ) 



PB = - 



1 0\ 

1 i 1 

4 \ Q 1 

1/ 



(34) 



in the {riB^riBi} = {0,0}, {1,0}, {0,1}, {1,1} basis. The 
von Neumann entropy of this is 



S(pb) = -tr(/5e log 2 Pb) 



(35) 



Therefore, according to Zanardi's site entropy entan- 
glement measure our example state contains two ebits 
of entanglement: one in the spin degree of freedom, the 
other in the spatial degree of freedom. This will be dis- 
cussed at length later in this paper. By contrast, as we 
have seen above, the Schlicmann measure gives zero en- 
tanglement for this state. We give in Appendix A an 
explicit construction showing that two qubits may be 
teleported using this state, further supporting the en- 
tanglement value given by the Zanardi measure. 



B. Behaviour under unitary transformations 

One-site two-particle (local) unitary transformations 
As before, we apply the infinitesimal one-site, two- 
particle unitary transformation (1 — ieH) with H = 
Un A ^n A ^. We obtain 



PB 




(36) 



1/ 



We now express this density operator as a density ma- 
trix in the binary {n A f, n^i, ^Ibj.} occupation num- 
ber basis, and reduce it for side B by tracing out states of 
side A using combinations of n A \ = 0, 1 and n A [ = 0, 1 
since the number of particles on site A is 0, 1 }> 1 }> or 2. 
Thus we perform 

PB= 

tia-\— 0,l,n_4|— 0,1 

(n Ah n Al \(nBT,n B i\p\n BV n Bl )\n Ah n Al ), 

(33) 



Hence, unlike the Schliemann measure, to first order in e 
the site entropy measure is invariant under one-site two- 
particle unitary transformations. This is the correct be- 
haviour for an entanglement measure: we cannot gener- 
ate entanglement through a purely local unitary trans- 
formation. 

Two-site one-particle (non-local) unitary transforma- 
tions Let us apply the transformation generated by ( |24] ) 

to our example state = -js{ c A -\ + c rt)^( c ai + 
(32) t 

c s|)|0)- Tracing out site A, we obtain the reduced den- 
sity matrix for site B, 



PB 



/I 



V. 



1 

-2iet 



-2iet 
1 



(37) 



1/ 



in the {hb^ubi} = {0, 0}, {1, 0}, {0, 1}, {1, 1} basis. To 
first order in e this is not equal to the untransformed ps- 
Therefore, two-site unitary transformations can generate 
entanglement in the site entropy picture, even if they only 
operate on one (delocalized) particle. This conclusion is 
as we would expect. 
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C. Site entropy measure applied to a completely 
general state 

Bosonic particles. Let us now apply the site entropy 
description of entanglement to completely general two- 
particle, two-site states. Since the case of bosonic parti- 
cles is the most general, we consider it first. The state 
can now be written in terms of the w-matrix as 



E 

.,b€{l,2,3,4} 



w ab bibi\o) 



(38) 



where 1, 2,3,4 = Af, Aj, -BT, B{ and w is now a symmet- 
ric coefficient matrix. 

Transforming this to the site-spin occupation number 
basis {n^f, Uai, n B[\ and tracing out the states of 
site A, we obtain a reduced density matrix for site B of 
block diagonal form, where each block corresponds to a 
particular occupancy (0,1, or 2 bosons) of that site. 



"bosonic 

Pb 



Pb,o 



Pba 



(39) 



PB,2, 



This is a 6 x 6 matrix, rather than the 4 x 4 pB we previ- 
ously obtained for the two-fermion state, because Bosc- 
Einstein spin statistics permit the extra site-i? double- 
occupancy states B1 £?| and B[ B{. 

The 'zero particles on site 23' component is 



Pb,q = \wu\ 2 + |w 22 | 2 + 4|wi 2 | 



(40) 



The 'one particle on site B' component in the B1, B{ 
basis is 

4|wi 3 | 2 + 4|w 23 | 2 4wi 3 u; 14 + Aw 23 w* 24 \ 
4^3^14+4^3^24 4|w 14 | 2 + 4|w 24 | 2 J [ ' 



Pb,i 



Finally, the 'two particles on site £?' component in the 
{B] B[,B] B],B[B[} basis is 

(4|uj 34 | 2 2w 34 w^3 2w 3i wl 4 \ 
2W34W33 |w 33 | 2 w 33 wl 4 (42) 
2wt, 4 W 44: W33W44 \W4 4 | 2 J 

Fermionic particles Obtaining an expression for ps 
for a completely general fermionic state is simply a mat- 
ter of applying the Pauli exclusion principle to p b ° sonlc . 
Under Fermi-Dirac statistics, the only possible two- 
particle state on site B is B f B {, meaning that the 
two-particle part of ps is reduced to the lxl submatrix 



(A 



'2. B J fori 



(4|W 3 4| 



(43) 



Similarly, the only possible two-particle state on site A 
is Af AJ,, meaning that the probability of zero particles 
on site B is given by 4|u;i 2 | 2 . Hence the zero-particle 
part of p h ° s °™ c is 

(/5 ,s)fcrmionic = (4|wi 2 | 2 ) (44) 

The one-particle part of p b ° somc is by definition not 
affected by Pauli exclusion, therefore 



(Pi,b) 



fermionic 



(Pi,b) 



bosonic 



(45) 



D. Relationship to Wootters tangle 

The origin of the Wootters entanglement measure is 
now readily understood. It is simply the von Neumann 
entropy of the one-particle part pf of the reduced density 
matrix in the occupation number representation for site 
B. Wootters's 'general state' equation (0), where the kets 
represent \<ja ob), can be rewritten in the occupation 
number basis \ua \ u aI «-bT ^b4) as 



a|1010> + 6|1001) + c|0110) + d|0101). 



(46) 



Tracing out site A yields the following reduced, cor- 
rectly normalized, density matrix for site B in the BT, B[ 
basis: 



|a| 2 + |c| 2 ab* +cd* 
PB ~ 1 a*b + c*d \b\ 2 + \d\ 2 



(47) 



with eigenvalues 
1 



1 



(1 - Vl-4|ad-H 2 ), g(l + Vl-4|ad-H 2 )- (48) 

Applying the simplifications r = 4\ad — bc\ 2 and x — 
2(1 + \/i — t) these reduce to 1 — x, x. Thus the entropy 
of pb is 

S( PB ) = ~{x\og 2 x + (1 - x)log 2 (l - a:)) (49) 

which is identical to the Wootters result for entanglement 
given in (^). 

V. SPIN-SPACE ENTANGLEMENT TRANSFER 
A. Omar et al. thought experiment 

Since we have argued that Zanardi's approach gives 
a correct view of entanglement in all circumstances, we 
can use it to analyze situations in which there is spatial, 
as well as spin, entanglement. A particularly interesting 
system of this type was introduced recently by Omar et 
al. [|l3| . They consider an apparatus which takes as its in- 
put two pairs of particles, A and B, each pair maximally 
entangled in some internal degree of freedom (e.g. spin), 
and transfers some of that entanglement to the spatial 
degrees of freedom of the particles. This is achieved by 
passing one particle from each pair through a beam split- 
ter on one side of the apparatus, and doing likewise with 
the remaining particles from each pair through another 
beam splitter on the other side of the apparatus (see Fig- 
ure ||). The two sides are labelled 1 and 2. Use of the 
site entropy measure enables us to understand better the 
process of entanglement transfer. 

Side 1 of the apparatus. First, let us consider the in- 
put state to the apparatus, and its entanglement accord- 
ing to the site entropy measure. This state is 



V=2 (a+A1 



± a 



All a A21; 



( a BlT a 



T"B 2 1 ^ a Bll a B 2 |) 



(.50) 



8 



Henceforth we will consider the case where all four par- 
ticles are fermions, and the above product state consists 
of triplets (described as the '++ case for fermions' in 
If we write this in the occupation number repre- 
sentation, and then trace out side 2 of the apparatus, we 
obtain the following reduced density matrix for side 1 of 
the apparatus: 



the apparatus: 



Pi. 



using the reduced basis for side 1 




(51) 



\n L i]n Ln n R1] n Rn ) = {|1100), |0011), 
|0110), |1001), |1010), |0101)}. 



(52) 



This state has two ebits of entanglement. Examining 
equation (|5l]), we see that this entanglement is carried 
entirely in the bottom right part of the density matrix, 
which corresponds to single-occupancy states which differ 
only by the spin. Therefore, this entanglement is purely 
spin entanglement. 

It is a straightforward exercise to show that the site en- 
tropy measure gives the same total entanglement between 
sides 1 and 2 of the apparatus (two ebits) for the input 
and output states. This must be so since the operation 
of each beamsplitter is local to its side of the apparatus. 
The unnormalized output state for 50/50 beam-splitters 
for our input state is given in E| as 

-i|L) 1 |L) 2 -i|i?) 1 | J R) 2 
-i(|£)i|i?} 2 + |i?)i|L) 2 ) 

+l{\An)x\m)2 + \mMnh) 

-i(|An) 1 |ATI)2 + |A|t)i|^IT) 2 ) 
+(|^tT)i|^4l> 2 + |A4|> 1 |ATT> 2 ) 



(53) 



where for example, |L)i indicates both fermions on side 
1 of the apparatus have passed into the left arm and thus 
necessarily have opposite spins, and |-ATI)i indicates that 
each particle on side 1 of the apparatus has passed into 
a different arm, with the particle occupying the left arm 
being spin up, the particle occupying the right arm being 
spin down. 

If we rewrite this in the occupation number represen- 
tation 



\nLi]n Ln n R i^n Rn ri L 2]n L 2iri R 2^ri R 2i) 



(54) 



Pi, out 
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* ?l 



(55) 



trace out side 2 of the apparatus, and renormalize, we 
obtain the following reduced density matrix for side 1 of 



which has entropy S(pi ]OU t) = 2, showing that the to- 
tal entanglement is unaffected by the operation of the 
apparatus. However, we can see from the fact that the 
double-occupancy top-left sector of this matrix is now 
non-zero that the system now contains spatial entangle- 
ment, because this state is now mixed in arm-occupancy 
number as well as spin. 

Single- occupancy and double- occupancy entanglements 
are additive. Since in ( |55"| ) there are no non-zero off- 
diagonal elements connecting the double-occupancy and 
single-occupancy sectors of the matrix, we can unam- 
bigously assign each eigenvalue to one sector, and hence 
divide the total entanglement into double-occupancy 
and single-occupancy parts. In this case, the double- 
occupancy sector has eigenvalues \ , and hence con- 
tributes 0.5 ebits to the entanglement. The single- 
occupancy sector has eigenvalues of j,0, j, \ and hence 
contributes 1.5 ebits. It is clear from these definitions 
that the single-occupancy and double-occupancy entan- 
glements will always sum to the total entanglement, pro- 
vided that the off-diagonal elements connecting the two 
sectors are zero. The single-occupancy entanglement is 
a form of spin entanglement, since the single-occupancy 
states do not differ in the spatial distribution of parti- 
cles between the arms. Likewise, the double-occupancy 
entanglement is a form of space entanglement, since the 
double-occupancy states do not differ in their m s values. 
However it is not obvious that these are the most general 
forms of spin and space entanglement, since for exam- 
ple, the double-occupancy entanglement does not take 
account of the spatial states in which each arm contains 
one particle. 

The distinction between spatial and double-occupancy 
entanglement is further clarified by the S x = spin mea- 
surements suggested by Omar et al. for their output 
state. They show that the spatial state produced by such 
a measurement (obtained with probability ^) involves a 
superposition of both double- and single-occupancy com- 
ponents. In this state, they show that the entanglement 
remaining between sides 1 and 2 is one ebit: since the 
spin state is now the same for all components and hence 
unentangled, this could be unambiguously described as 
spatial entanglement. 

Left arm of side 1 of the apparatus. It is instructive 
now to reduce further the input and output density ma- 
trices to those for just the left arm of side 1 of the appa- 
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ratus. For the input state, this is 



V 



in the basis 



\nmn Ln ) = {|11),|00),|01),|10)} 



(56) 



(57) 



which has entropy S^/SiL.m) = 1- In the same basis, the 
reduced density matrix for the output state is 



PlL.out 



n 



v 



v 



(58) 



which has entropy S (/Six, out) = 1-81, showing that the 
action of the beamsplitter on side 1 of the apparatus 
has introduced an additional 0.81 ebits of entanglement 
between the left arm of side 1 and the rest of the system, 
in addition to the 1 ebit of entanglement already present 
between those two subsystems. 

Operator-sum representation for spin-space entangle- 
ment transfer. It is possible to find an operator-sum 
representation for the spin-space entanglement transfer 
within the left arm of side 1 of the apparatus that we 
have discussed above. An easy way to do this is to make 
the following isomorphism between the spin states of two 
qubits A and B, and the occupation numbers for the spin- 
up and spin-down single particle states of the left arm: 

Wactb}= {|TT), IU), lit), III)} 

{n 1L1 n 1Ll }= {|00),|01),|10),|11)}. (59) 

We then find that the action of the Omar interferometer 
in transforming pi^^n to /5iL, ut can be represented by 
the action of the depolarizing channel Q on piL.m with 
probability p = | . 



B. Division of entanglement into single- and 
double-occupancy parts 

When is it possible to divide entanglement unambigu- 
ously into single- and double-occupancy parts? As we 
can see from the above treatment of the Omar appara- 
tus, it is when the reduced density matrix for the subsys- 
tem whose entanglement we are considering has sectors 
corresponding to single- and double-occupancy, with no 
off-diagonal elements connecting them. Such a division 
is possible whenever the total system contains a definite 
number of particles: there are then no off-diagonal den- 
sity matrix elements connecting states of the subsystem 
having different numbers of particles. This is the rea- 
son why there are no elements connecting the different 
sectors of pn, in equations (B6J) and (J58|). 



We note with emphasis that the situation for the re- 
duced density matrix for side 1 of the Omar apparatus 
in equation (|5^) is fundamentally different. All the ba- 
sis states contain the same total number of particles. Its 
block-diagonal form is due to a combination of factors: 
the spin symmetry of the system (which causes those 
elements connecting states on side 1 and 2 with a dif- 
ferent total m s value to be zero), and the use of 50/50 
beam splitters, which prevents any products of the form 
|L)i|j4) 2 from appearing in the output state. 



APPENDIX A: TELEPORTING TWO QUBITS 
USING AN EXAMPLE DELOCALIZED STATE 

Protocol design. Consider again the delocalized state 
in equation (|20j). Since the Zanardi measure says it con- 
tains two ebits of entanglement, we should be able to tele- 
port two qubits of quantum information using it. Clearly, 
since the two ebits are spread across spin and spatial de- 
grees of freedom, we shall need to modify the original 
protocol somewhat. How could we do this? 

Switching into the anthropocentric language of 'Alice' 
and 'Bob', a concise description of the protocol for tele- 
porting one qubit described in Q is as follows. We sep- 
arate the two subsystems of our entangled system which 
will act as the channel for quantum information, giving 
one to Alice and one to Bob. We then perform a CNOT 
on the qubit whose state we wish to teleport (the 'source 
qubit'), and Alice's system, using the source qubit as the 
control line. We then perform a Hadamard transform on 
the source qubit. Alice's qubits are now in a superpo- 
sition of states, each of which corresonds to the target 
qubit being in the same state as the original state of the 
source qubit up to a unitary transform. Alice performs 
a measurement of the state of her two qubits, thereby 
projecting the target qubit into one of these states. The 
protocol is completed by Alice sending Bob two bits of 
classical information describing which measurement re- 
sult she obtained, enabling him to rotate the target qubit 
into the correct state. 

The key to teleporting via the delocalized state (^) lies 
in recognizing that the two ebits in the delocalized state 
are equivalent to two pairs of qubits, each of which is 
maximally spin-entangled ('channel pairs'), and making 
the following isomorphism: 

{o A a B }= {|TT), ITI), lit), III)} 

{n A] n Ai }= {|00),|11),|10),|01)} (Al) 

This connects the occupation numbers of the single par- 
ticle states of Alice's site to the states of Alice's channel- 
pair qubits in the spin-only representation. 
Recall that a CNOT performs 

ITT) - IT!), IT!) - ITT), III) - IIT), III) - III) (A2) 

i.e. we flip the second qubit in a basis state iff the state of 
the first (control) qubit in that basis state is 'up'. What 
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does a CNOT on one of Alice's two channel-pair qubits 
look like after applying the above isomorphism? Using 
this basis for the states of one of the source qubits (C) 
and Alice's site (A): 



{\ncincin A1 n Al }} 



{|1000), |1010), |1001), |1011), 
|0100), |0110), |0101), |0111)} 

(A3) 



we obtain the following unitary transformation for the 
first 'virtual' qubit: 

/0 1 0\ 
1 
1 



u 



00100000 
~ 00000000 
00000000 
00000000 
\0 0/ 



The action of this is thus: 

|1000) <-> 1 1010), |1001) <-> |1011), 
|01nAt^Aj.) unchanged. 



(A4) 



(A5) 



Referring to the isomorphism in (|Al| ) we see that this U 
flips the first 'virtual' qubit in Alice's half of the delo- 
calized state iff the control qubit C is spin-up. Similar 
considerations lead to a similar unitary transformation 
for the second 'virtual' qubit. We also note that since we 
are teleporting two qubits, we need to send four classical 
bits to complete the protocol. 

Protocol implementation. The two CNOTs described 
above will clearly allow us to exploit the two ebits of 
entanglement present in the delocalized state. However 
some consideration needs to be paid to how we can imple- 
ment these CNOTs. Considering again the first 'virtual' 
qubit, the Hamiltonian we can use to generate equation 
(A§) is 



qubit = |10)cc(10| I |00)aa(10 



|10)aa(00| + |11)aa(01| + |01)aa<H| 



-K.c + 1) c AT + c A] 



(A6) 



In the first expression we have used the bases 
\ncfnci) and \n A ^n Ai ). 



(A7) 



In the second expression we have rexpressed the projec- 
tors for the occupation number state of site A in second- 
quantized notation. Similarly, the Hamiltonian generat- 
ing a CNOT on the second 'virtual' qubit is 



1 



1) 



c \] c \[ + C M C A] )-(A8) 



Hind qubit 

Neither of these Hamiltonians conserves particle num- 
ber, thus we need to introduce a coherent source/sink of 
particles to the system. We shall see in a moment that 
we can easily do this for bosons. Introducing a system D 
which acts as a particle source/sink, Hi st qu t>it becomes 



Hut qubit = |10)OC<10|( C^Cu + cJjCAt )• (A9) 

At this point we face a problem. By changing the 
number of particles in system D as a consequence of our 
CNOT, we are introducing new correlations between the 
states of subsystems A and D. This is thus a type of de- 
coherence affecting the entanglement of the 'carrier-pair' 
AB. This is clearly unavoidable in a real- world system, 
but we can show that for bosons, by choosing a suitable 
initial state for subsystem D we can minimize this deco- 
herence to a negligible level. We seek to put system D 
in an approximate eigenstate of the creation and annihi- 
lation operators, so that they leave it unchanged and no 
decoherence of the entanglement in the AB carrier pair 
occurs. A suitable choice is the coherent state 



a) D 



exp(-i|a| 2 



E 



(n!)V2 



\n) D . 



(A10) 



It is well known that this state is an eigenstate of the an- 
nihilation operator, a fact which suits our requirements 
perfectly, but it is not an eigenstate of the creation op- 
erator. However, as the mean number of particles \a\ 2 
in the coherent state asymptotically approaches oo, the 
state asymptotically approaches an eigenstate of the cre- 
ation operator. It is important to note that this method 
for coherently producing a non-number-conserving inter- 
action applies to bosons only. For fermions, Pauli exclu- 
sion prevents us using such a simple approach and there 
is no analogue of the coherent state available within the 
Hilbert space. 
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FIG. 1: Schliemann entanglement r\ of all Bell states vs overlap S = \{cf> a \<j>b) \ of single particle states. 
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FIG. 2: The spin-space entanglement apparatus used in the Omar et al. thought experiment (reproduced from [[l3[) 
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